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1. INTRODUCTION 

The methods of feedback control are widely used 
in the modern physics, but still they are not very 
popular in quantum optics. Very often this "cy- 
bernetical" approach does not demand involve- 
ment in very complicated physical devices and can 
be arranged in a trivial nonlinear system [1] . 

We apply this technique to control the energy of a 
two- level atom in the optical external field E(i) in 
the frame of the so-called "semi classical model" of 
the atom-field interaction that describes a single 
quantum two-level atomic system (all other levels 
are neglected) with classical electromagnetic field. 

Recently other authors studied the control of two- 
level atoms in the frame of open loop-ideology 
when the controlling field was known a priori. It 
allowed obtaining the different forms of atomic 
energy spectra, producing ir- and ~ pulses [2], 
taking special non-constant shapes of external 
field [3] etc. 

The main feature of the model proposed in this 
article is that, it is based on the closed- loop 
approach. This means that we do not initially 
define the dependency of the field on time, but 



restore this function for every moment from the 
current values of the probability amplitudes of the 
atomic ground and excited levels. 

The closed-loop (feedback) scheme for the inter- 
action of two-level atom with external field can 
be realized in different models. The most famous 
and fully developed is the approach based on 
master equations in its both main variants: the 
Markovian feedback model [4] and the so called 
Bayesian feedback [5] (the later model was pro- 
posed by Wiseman in his comparative analysis of 
both these models in [6]). The Bayesian ideology 
is more closely related to our approach because 
this closed-loop control is constructed directly on 
the estimation of the system state. 

Another approach is to construct the control 
scheme for a single atom for the quantum control 
field [7]. In this paper we discuss the classical 
control field and we do not apply special restric- 
tions on its shape i.e, the optical field shouldn't 
be sinusoidal as in [8] or have other special time 
dependency. Thus, our scheme of classical feed- 
back proposed here is similar to the traditional 
variant of control theory in the form of speed- 
gradient (SG) method [9], when input variables 
change proportionally to the speed-gradient of 



appropriate goal function. We use the standard 
notations following [10], but in our model the 
optical field plays the role of a control signal u(t) 
for closed-loop or feedback control scheme. 

For this purpose we use the real positive goal 
function Q, measuring how far at the moment we 
are from the desired state of the atom. As a result, 
we calculate the control signal u(t), i.e. we restore 
the shape of the electromagnetic filed E(t) to keep 
the atom at the upper level. 

In the second section of this work, we construct 
the feedback control model for the single two-level 
atom in external controlling optical field. Then, 
in the third section, we apply feedback speed 
gradient scheme to the non-decay case. 

2. TWO-LEVEL ATOM IN CONTROL 
OPTICAL FIELD 



Let's consider the interaction of an optical field 
E(t) linearly polarized along the x-axis with a 
two-level atom. 




where p ab — p ba — e(a\x\b) is the matrix element 
of the electric dipole moment. We neglect the 
decay of the levels. We express the electric field 
as 

E(t) = E Q u(t), 

where Eq is the amplitude and u(t) is the dimen- 
sionless control signal. The equations of motion 
for the amplitudes C a and C b may be written as 

C a = ~Luo a C a + if2_Ru(t)e _t0 Cb; 

C b = -mj b C b + tQ R u(t)e^C a , 

where the "Rabi frequency" is defined &s £l R = 
\ pba ^ E ° ; and <f> is the phase of the dipole matrix 
element p ba = |pb a |e^. 

To solve for C a and C b , we write the equations of 
motion for the slowly varying amplitudes as: 

c a = C a e ,aJa * ; c b = C b e LUJbt , 

then 

c a = Ltl R u(t)e-^c b e' ut ; 

where u> = oj a — tu b is the atomic transition 
frequency. The phase <fr can be excluded from the 
system, if we put c b = c b e~^: 

cl - Ln R u(t)e- LWt c a 

Later for simplicity we will denote c b with c b , then 
finally: 

6a = Ln R u(t)e LWt c b (1) 
ct = Ln R u(t)e~ lwt c a (2) 



Fig. 1 . Interaction of a single two-level atom with 
an optical field. 

Let | a) and \b) represent the upper and lower level 
states of the atom, i.e. they are eigenstates of the 
unperturbed part of the Hamiltonian Hq with the 
eigenvalues: Ho\a) = hu a \a) and Ho\b) = hu> b \b). 
The wave function of a two-level atom can be 
written in the form 

\iP(t))=C a (t)\a)+C b (t)\b), 

where C a and C b are the probability amplitudes of 
finding the atom in states \a) and \b), respectively. 
The corresponding Schrodinger equation is: 

with H = Ho+Hi, where Ho and Hi represent the 
unperturbed and interaction parts of the Hamil- 
tonian, respectively [10]: 

Ho = hw a \a)(a\ + huj b \b){b\ ; 

Hi = -(p ab \a){b\ + p ba \b){a\ s jE{t), 



Now let's suppose that we have the initial condi- 
tions: 

c Q (0) = ; c 6 (0) = l (3) 

and our goal is to stabilize the atomic system at 
the upper level: |c a | 2 = 1. 

3. SPEED GRADIENT METHOD FOR 
PROBABILITY AMPLITUDES CONTROL 

We have not yet specified the time-dependent 
function u(t). To find it, we apply the speed 
gradient (SG) method [9] to control the system 
behavior. 

In this approach, the control action is chosen in 
the maximum descent direction for a scalar goal 
function. 

The goal in the control process is a smooth scalar 
function Q with the limit relation 

lim Q(x{t),t) -> 0. 



The purpose of the SG method is to minimize the 
goal function 



\{\c a \ 



1 



(4) 



where |c a | 



SG represents the control signal u with the time 
derivative of the goal function Q. 

The underlying idea of SG method is that moving 
along the anti-gradient of the speed Q provides 
decreasing of the goal function. In our case the 
control signal space is 1-dimcnsional, thus we 
reduce our gradient to the partial derivative with 
respect to u. In the case of proportional feedback 
with some positive coefficient T, it is defined in 
the form: 



u = -r 



dQ 
du 



(•5) 



Thus 



«(t) = Lm R (\c a \ 2 - 1)( 



Milt * 

e e a e b 



vt c b <)6) 



Putting value of u(t) from Eq.(6) in Eqs.(l) and 
(2), we have the following system of equations: 



Thus, we can simplify the system (7), putting 

Pb = 1 — Pa- 

The system (7) has two equilibrium (fixed) points: 
(p a ,p+,P-) = (0,0,0), (1,0,0) 

On the Figs. 2,3 we demonstrate the result of 
our control procedure for: T = 0.1 sec, fl R = 
10 2 sec -1 and uj = 10 3 seer 1 . 

On Fig. 2 we show the solution of Eq.(7a). 
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Fig. 2. The density matrix element p a (t) for the 
control procedure (4)- (5) 



liwt * 2 I 1 2 

e c a c b ~ c a|Q> 



c a = Ttt 2 R (\c a \ 2 
c b = Ynl{\e a \ 2 - l) (c b |c a | 2 - eyelet) 
Now suppose that 

Pa = C a C* a — \c a \ 2 ; Pb = C b cl = \c b \ 2 ] 

up— — e c a c b — e c a c b , 
p+ = e c a c b + e c a c b . 

Hence we have the following four equations: 



p a = 2r^ 2 i ( Pa -i)[( 



2 1 

'pi. -pi 



PaPb 



p b = 2vn' R \ Pa -i 



PaPb 



P\-Pi 



P+ = up- ; (7) 

P- = -2™ 2 R (p b - p^j ( Pa - - LU P+ . 

Also from Eq.(6) the control signal u(t) becomes 

u (i) = -rr> fl (p a -i)p_ (8) 

With initial conditions p a (0) = 0, p b (0) = 1 we 
have 

Pa + Pb = , 

that means in fact: 

|c Q (t)| 2 + |c fc (i)| 2 = l, 

which is the simple statement that the probability 
to find the atom in one of its states \a) or \b) is 1. 



On Fig. 3 we show the control signal(8). 
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Fig. 3. The control signal u(t) for the system (7). 



4. CONCLUSION 

The SG algorithm can be easy applied to establish 
feedback control for the probability amplitudes of 
two-level atom. 

This scheme can be modified if we take into con- 
sideration the decay of the atom levels, because 

in this case the goal Q = ^|c a | 2 — ij is not 
achievable for SG algorithm in principle. For this 
purpose we will redefine the goal function Q. 
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